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ABSTRACT 

A basic method is presented for the estimate of the overall mechanical 
response of solids which contain periodically distributed defects (inhomogeneities, 
regions undergoing inelastic flow, voids, cracks, etc.)* This method is then 
applied to estimate the shape and growth pattern of voids that are periodically 
distributed over the grain boundaries in a viscous matrix. The interaction effects 
are fully accounted for, and the results are compared with calculations for a 
single void in an infinitely extended viscous solid, by Budiansky, Hutchinson, and 
Slutsky. Then, for a polycrystalline solid that undergoes relaxation by grain 
boundary sliding, the relaxed moduli are obtained, again fully accounting for the 
interaction effects. Finally, the overall inelastic nonlinear response at elevated 
temperatures is discussed in terms of a model which considers nonlinear power law 
creep within the grains, and linear viscous flow in the grain boundaries. 
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1. INTRODUCTION 


The Inelastic response of polycrystalline solids stems from a variety of micro- 
structural changes, depending on the temperature regime, as well as the stress 
history. At temperatures sufficiently below 50% melting point, the rate effects are 
not dominant. The deformation consists of plastic slip on crystallographic planes, 
accompanied by the accommodating elastic lattice distortion. At higher temperature 
regimes the rate effects become significant, and the intracrystalline flow can be 
modelled adequately by a power law. At higher temperatures, creep effects are the 
major components in the overall response. In this case, adequate micromodelling 
involves a power law flow within the grains, accompanied by linearly viscous grain 
boimdary sliding. 

Various micromechanical defects that may arise in the course of deformation, 
contribute differently to the failure mechanisms during different temperature 
regimes. At low temperatures, voids are generated because of plastic flow at second 
phase particles, and this may lead to a reduction in ductility. At higher tempera- 
tures, on the other hand, voids are nucleated on grain boundaries, and grow in re- 
sponse to the applied load, as the solid creeps. Depending on the load level and 
the temperature regime, the mechanism of such void growth varies. For example, 
cavity growth is essentially crack-like, when surface diffusion is much slower than 
the grain boundary diffusion, whereas at a high stress level the cavity grows essen- 
tially by intragranular power law diffusion. 

Under NASA-Lewis sponsorship, theoretical and experimental work has been initi- 
ated at Northwestern University on the micromechanical modelling of nonlinear con- 
stitutive relations of superalloys at various temperature regimes, addressing all 
the above-mentioned microscopic features. The present report summarizes some of the 
theoretical results on the growth of voids in viscous metals, the effects of grain 
boundary defects on the overall response of the polycrystal, and, finally, the over- 


114 



all creep response of the polycrystal. Elasto-plastic (rate-independent) modelling 
is discussed in a separate report; Iwakuma and Nemat-Nasser (1982). 

The calculation of the overall response of the polycrystal is based on some 
fundamental results on the effect of periodically distributed defects (inhomogenei- 
ties, regions undergoing inelastic deformation, etc.) on the overall response of the 
solid; Nemat-Nasser et al. (1982). These results are first briefly discussed, and 
then applied to the estimation of the shape and growth pattern of voids that are 
periodically distributed over the grain boundaries in a viscous matrix. The inter- 
action effects are fully accounted for, and the results are compared with calcula- 
tions for a single void in an infinitely extended viscous solid, by Budiansky, 
Hutchinson, and Slutsky (1982). Then, for a polycrystalline solid that undergoes 
relaxation by grain boundary sliding, the relaxed moduli are obtained, again fully 
accounting for the interaction effects. Finally, the overall Inelastic nonlinear 
response at elevated temperatures is discussed in terms of a model which considers 
nonlinear power law creep within the grains. 
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2. FORMULATION OF THE BASIC PROBLEM 

Consider a solid containing periodically distributed sets of Inhomogeneltles 
such that It can be regarded as a collection of Identical unit cells. Let D be a 
typical cell of volume V and exterior surface S. For simplicity assume that D Is a 
parallelepiped of dimensions measured along the rectangular Cartesian coordinate 
axes x^, 1 1,2,3. The results also apply to a single cell subjected on Its bound- 

ary to suitable displacement or velocity fields. 

Neither the matrix nor the Inhomogeneltles are required to be linearly elastic 
or rate-independent, but, for the intended applications, only small strains and 
rotations are considered. 


To be specific, let g be the Cauchy stress and set 

■“"IJ ■ 

where repeated indices are summed over 1, 2, 3, dg Is the stress Increment, and D *= 
D(a) is the instantaneous compliance which may or may not depend on stress. For 
rate-independent applications, e in (2.1) is the strain tensor. For rate-dependent 
cases, on the other hand, e Is the strain rate tensor. For example, for non-linear 
creep, the strain rate is e = F(g). In this case we consider the Incremental rela- 


tion 




In particular, if power law creep is assumed, e' •= nJ ct', we obtain 
-nr. . , "'^11‘^kAi , , 

where prime denotes the devlatoric part, and 


( 2 .^ 

is the effective stress; in (2.3) n Is a positive number and q is a dimensional 
parameter. If an inhomogeneity is linearly elastic or linearly viscous, then D in 
(2.1) would be a constant tensor with suitable usual symmetries. 
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Let C be the Inverse of D and rewrite (2.1) as 




(2.5) 


Again, C may be a function of q. 

Assume now that the displacement (velocity) field u° is prescribed on S in 
such a manner that the average strain (strain rate) field is obtained. Let the 
corresponding average stress field be q. Consider an incremental change, du*^, in 
u°, which produces the increments de° and dg in the average strain (strain rate) 
and stress fields, respectively. We seek to calculate the overall moduli C*, defined 
by 


where, in general, C* depends on the average stress a, as well as on the micro- 
structure. 

Within the unit cell, neither the stress increment nor the strain (strain rate) 
increment is uniform. Let there be M inhomogeneities, r = 1,2,...,M, and set 






+ dCkj(j)] 


in D - n 

r 


(2.7) 


where de is the perturbation strain (strain rate) field due to inhomogeneities; C is 
the modulus tensor of the matrix; and is that of the r^^ inhomogeneity. 

As has been shown by Eshelby (1957) for an ellipsoidal inhomogeneity in a 
linearly elastic, unbounded solid, the nonhomo geneous body may be replaced by a 
homogeneous one, provided that suitable transformation strains are prescribed in the 
ellipsoid. In this case, the transformation strain tensor is constant. For period- 
ically distributed inhomogeneities, or when the inhomogeneity is not ellipsoidal, 
the transformation strain tensor is no longer constant. The basic concept, however, 
still applies, and can be quite effective, as shown by Nemat-Nasser and Taya (1981) 
and Nemat-Nasser et al. (1982). 
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Hence, in place of (2.7), one writes 




( 2 . 8 ) 


where dg*^ is zero in and seeks to express this transformation strain (strain 

rate) increment in terms of de^. This is done by the use of the Fourier series 
representation of the incremental fields, as has been discussed by Nemat-Nasser et 
al. (1982). The final results for the present case are as follows: 


a5 “ [c - c 

Tc£mn k£pq k£pq pqmn’ 

I / 8cjJ(x')eM-<5-?'>dx-. 

* n =0 r=l n 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


(no sum on j ) , i = 


( 2 . 12 ) 


and where k,£,m,n,j = 1,2,3. In (2.11), the fourth order tensor (£) depends on 

the matrix modulus tensor C. For an isotropic matrix. 


^ijk£ “ ^^'^ik*^j£ ^i£^jk^ 


(2.13) 


®ijk£^l^ “ ^j^‘^i£^k‘‘’'^ik^£^ ^i^'^j£^k'*‘ ^jk^£^ " 1- v ^i^j^k^£ 1- v ^i^j'^k£’ 


(2.14) 


5i -= 5^/5. 


^ = <^k^k> • 


In (2.13) and (2.14), X, y, and v are material parameters for the matrix, which may 
depend on stress g. For a linearly elastic matrix, these are the usual Lamd con- 
stants and Poisson's ratio, respectively. In the general formulation that will fol- 
low, we shall assume an anisotropic matrix . In Section 3, however, we assume an 
isotropic matrix, and hence use (2.14). In Section 4, on the other hand, a non- 
linear creep law is considered, and this makes the tensor C dependent on the current 
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stress state; then (2.14) cannot be used, and hence a more general expression is ob- 
tained . 

Let be the volume fraction of the r^^ inhomogeneity, 

= V^’ ^2.15) 

xir 

where is the volume of and denote by de the average value of de* taken over 




dg*^ = ^ / de*^(x)dx; 


(2.16) 


r n. 


note that de*^ is zero outside of Averaging (2.8) over D, and using (2.6) we_ 

obtain 


M 




(2.17) 


where the notation 




(2.18) 


is used. We now substitute from (2.11) into (2.9), average the resulting equation 


over to arrive at 


where 


'jk jkmn’^mn ^ ®jkmn 


r n 


s=i n 


s 


( 2 . 20 ) 


It has been shown by Nemat-Nasser and Taya (1981) that good accuracy is obtained 
if the transformation strain (strain rate) increment in the integrand in (2.19) is 
replaced by its average value. This then leads to 


M 


^r ‘^^jk ^jkmn^mn " ^jkmn^mn* 


( 2 . 21 ) 


where 
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( 2 . 22 ) 


“P 

Equations (2.21) are now solved for B , results substituted into (2.17), and since 

mn 

de° is arbitrary, the following general result is obtained: 


ijkJl " ^ijki “ ^ijinn ^s^ “ ^mnkil ^r^ ’ 


(2.23) 


s=l r=l 


where 6 is the Kronecker delta. 

In (2.23) the tensors C, A^, and S""* may depend on the stress, a, in the matrix 
as well as in each corresponding inhomogeneity. The estimate of the stress variation 
throughout the solid is indeed a formidable task. For our purposes, it seems ade- 
quate to use the overall average stress a instead. Then the overall stress-strain 
(strain rate) relation can be obtained incrementally with the aid of (2.23) and 
(2.6). Some specific results are presented in subsequent sections. On the other 
hand, when necessary the local strain (strain rate) increment in, say, can be 
obtained from (2.9), 


and hence the local stress increment can be estimated from (2.7), 


daj.(x) = cj., , Af. de*’^(x) in . 
ii ~ iik£ TcJlmn mn - r 


In a similar manner, the stress increment within the matrix can be obtained from 
(2.11) and (2.7). 


120 



3. GROWTH OF PERIODICALLY DISTRIBUTED VOIDS IN VISCOUS METALS 


3.1 Background 

At elevated temperatures, voids are nucleated at grain boundaries in polycrys- 
talline solids. Depending on the deformation and temperature histories, the arrange- 
ment of these voids relative to the orientation of the principal applied stresses 
can vary considerably. For example, experiments show that voids can be concentrated 
on grain boundaries perpendicular to the direction of maximum tension, see, e.g., 
Garofalo (1965). For superalloys that are plastically deformed at room temperature, 
on the other hand, Dyson et al. (1976) have shown and Kikuchi and Weertman (1980) and 
Saegusa et al. (1980) have conclusively verified that after annealing, voids are 
generated at grain boundaries parallel to the direction of maximum tension. The 
mechanisms giving rise to the formation of these cavities are different, but their 
presence has similar adverse effects on the life expectancy of the structural com- 
ponents. An account of diffusive cavitation in polycrystalline solids is given by 
Chuang et al . (1979) and by Argon al . (1981); see also Rice (1981). Here, how- 
ever, a different approach is used, which considers the growth of periodically dis- 
tributed cavities within a viscous metal. We make contact with the work by 
Budiansky ^ al . ( 1982) who examine the growth of a single cavity in an unbounded 
viscous medium, as well as with an earlier contribution by McClintock (1968) on the 
same subject. 


3.2 Formulation 

For a linearly viscous matrix , we have 
^Ij “ ^ijki^kJl’ 

where 


^ijki, “ ^^^^ijkZ l-2v '^lj\il^’ 




ijk£ ijk£ 2y' ijkil 1+ v ij kil 


(3.1) 


(3.2) 

(3.3) 
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In (3.2), is the fourth order identity tensor; in (3.1) e is the strain rate . 

Consider a unit cell of dimensions A^, and let it include an ellipsoidal void of 
(principal) semi-axes a^, oriented along the coordinate axes x^, i = 1,2,3. Define 


. B = — , Y = 7-. f = 


4tt 


a^ ’ ^ Aj ’ ^ Aj ’ 3 AiA2A3* 


(3.4) 


Since (3.1) is linear, all the Incremental relations in Section 2 can be re- 
placed by the total ones, i.e. all the relations apply if de'^, da,..., are replaced 
by £°, a,..., respectively. 

From (2.8) it follows that^ 

e* = e^ + e (3.5) 

within the void, and from (2.21) we obtain 


^ij “ ^^IjkJl " 

where, in view of (2.22) and (2.14), 
c = _ i o ^ ~ g 
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"T 

1- V 
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^2211 
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.1- V 

1 

s.. 

4- 

V 

s„. 

s _ - - “ 
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s„ 
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1- V 

7 


1- V 

2’ 
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2 '■2 

3 


1 - \; 


s = — — - — s + — — s 

^3322 1 - V 7 ^ 1 - V 3 ’ 


r. s — — <> + 

‘’3311 1 -V 8 ^ 1 - V 3 * 


i«3 + Sj) - 


^211 - - ihh- 


The infinite series S^ = S^(a,B,Y>C>f ) (3.7) is defined by 

+00 

S. “ ^ P(n)h.(5)> 1»2,...,9, 


n =0 
P 


(3.6) 


(3.7) 


(3.8) 


where 


"since M = 1, the superscript r = 1 is dropped. 
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(3.9) 


hg = a,)\ « a^K,)^, hg - (535^)2, hg = aj^y. 


and, for an ellipsoidal inhomogeneity , 


p(n) = f 


9 (sin n ~ Ti cos T))^ 


n ?* 0, 


r on 2 gn- ^ 1*5 

n - 27rR[n| + (— ) + (— ) ] . 


r 2Trf 2 ^ 


In view of (3.5), the shape change can be defined by 

• • • • 

a, a V„ 

—1 = Z* a Z* ~ « E* — = 7* 

a, ^ 11 ’ a, ^ 22 ’ a, ^33’ V„ ’ 

1 2 3 n 

and we also note, from (3.4), that 


is. a ”1* Z* 
a 22 ^11 ’ 


*® e* “ E* 

B 33 11 ’ 


r -0 - ^0 


'22 '"11 


e w 

1 1 > 


i a pO _ pO 1 = 7* - pO 

c ^33 ^11 ’ f ^k 


To obtain the current dimensions and other geometric variables, we integrate (3.12) 
and (3.13) with respect to time. This, for example, yields 




l^'o 0 “00 


^22 " 


where the subscript 0 denotes the initial value. 

Since the transformation strain rate tensor, g*, characterizes the rate of 
change of the void geometry in accordance with (3.12) and (3.13), Eq. (3.6) relates 
the void change parameters to the overall strain rate tensor e°. To make contact 
with results of Budiansky et al. (1982), we relate the overall strain rates to the 
average stresses by 


^4.. = C* , „ e°„. 
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and note that, unlike the case of a single void in an infinitely extended solid con- 
sidered by Budiansky et al. (1982), here C* does not equal the matrix modulus tensor 
C. The overall modulus in the present case is obtained by specializing (2.23) or, 
equivalently, by equating the overall rate of energy loss per unit volume with the 


average rate of loss. This results in 


r* -OgO pOpO pOp* 

^ijki^ij^kJl ^ijki ij kZ ijkZ^ij^kZ 


(3.16) 


which, for M “ 1 and because of (3.6), implies (2.23). Since is arbitrary, 
(3.16) and (3.15) yield 


^ij ' ^ijkZ^^U ■ ^ ^kZ^' 

In the present case C is isotropic, Eq. (3.2), and if we introduce 


cO 




and eliminate between (3.17) and (3.6), we obtain 


(3.17) 


(3.18) 


1 +v 


1 -V -V 

-V 1 -V 

-V -V 1 
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11 


^ ^ “ ^1111 

"^1122 

“^1133 

0 
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► = 

“^2211 

1 - f - S2222 

“^2233 

0 

3 3 >' 


- “^3311 

“^3322 

1 - f - S 


3333 


( 7 * ) 
11 


< e* 

‘ 22f ’ 


7* 
'>33'' 


'12 


'23 


(1 - f - 


(1 - f - 2S2323 )e23» 


cO 

^31 


:* 

'z: 


(1 - f - 2S3j3^e2j. 

From (3.19) it follows that 
1- 2v 


(S»J + S»^ H- SO3) - {1 - - 833^3 - 833^^ - f}S *3 


+ {1 - s 


2222 ^1122 “ ^3322 


-f}< 


“ ^3333 ~ ^1153 “ 


1133 


2233 


22 

T* 

33 ' 


(3.19) 


(3.20) 


^The calculation is essentially the same as in Nemat-Nasser et al . (1982). 
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For the incompressible matrix, v = 1/2, and for (3.19) and (3.20) to yield non- 
trivial results, the matrix in the right-hand side of the matrix expression in (3.19) 
and the coefficients in the right-hand side of Eq. (3.20) should vanish.^ With 
V ■ 1 / 2 , this leads to 


1 - 3 S^ - - S3 + 2(S^^ + Sg + Sg) = f, 

1 - - 3S2 - S3 + 2(85 + + Sg) = f, ( 3 . 21 ) 

1 - Sj - S^ - 3S3 + 2 (Sg + S7 + Sg) = f. 


and, if only the infinite series S^ , S 2 , and Sg are retained, from (3.7) we deduce 
that 


nil 

- s 

- s 

r 1 - 2v 

^2211 

3311 

“ 2(1- v) 

2222 

^1122 

^3322 

^ 1 - 2v 

“ 2(l-v) 

3333 

- s 

- s 

r 1 - 2v 

1133 

2233 

" 2(1- v) 


- £ - S, 

{1 - £ + Sj 

{1 - £ + Sj 


+ Sj + S 3 }. 

-S3+S3I. 

+S 3 -S 3 }. 


(3.22) 


With these and with v = 1/2, (3.19) yields 


er. 


■ 1 s° - i s° - i sO ) 

11 


3 11 3 22 3 33 

-22 

. " tly]. 

_ i cO + 2 gO _ L cO 

3^ir 3 ^22 3 ^33 

>^33. 


1 1^11 3 “22^ ^ ^337 


(3.23) 


where is the inverse of the matrix 


-Sun 

"^1122 

"^1133 



“^2211 

1 - f - S2222 

“^2233 

• 

(3.24) 


1 — f — Sj + S 2 + Sg 


1 - f + Sj - S^ + Sg 


1 - f + Sj +S 2 - S 3 


Equations (3.23) relate the void growth parameters to the overall stress components. 
In terms of the stress ratios 


'Numerical tests for spherical, cylindrical, and ellipsoidal geometries show 
that to within the accuracy of the estimate of the infinite series, these conditions 
are almost satisfied. 
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> 


(3.25) 


^22 *^22 


S° 


33 33 


'11 011 


'11 011 


one obtains 


3e* 

•^^11 

^11 

^^22 

S^’ — 
^11 

3e* 
^^33 


(2Tii-T,2 + 2Ti3) + (-T^,+2T,2 + 2Tj 3)<J> + ("Tn - 2Ti3>^, 


(2T21 -T 22 + 2 T 23 ) + (-T,i +2T„ +2T„)(J + (-T, , - T„ + 2T, . 


"21 ■ "^-^22 ^ ‘"••^23 


21 22 23 ' 


S^= (213,-132+2133) + (-T3, +2132+2X33)^ + (-'^ 31 -^ 32 ‘"^^ 33 >'^- 


11 


(3.26) 


Finally, the components of are obtained from (3.6) and (3.26). 


3.3 Numerical Results 

Table 1 lists the initial and the loading conditions for eight different cases 
which are reported here for illustration. It should be noted that even in high 
strength metals which undergo very small overall deformations, the local deformations 
close to inhomogeneities or at the tip of cracks can be quite large. For this reason 
in Fig. 1, results for rather large strains are included. This figure shows the void 
volume change as a function of the overall deformation measure, L/Lq or L^/L, for the 
indicated cases associated with Table 1. For comparison, an asymptotic and addition- 
al results of Budiansky (1982) are also shown. [These are read off the fig- 

ures in the published paper. In the final version of the present report, these will 
be recalculated in order to obtain a more accurate estimate of the effect of perio- 
dicity as compared with a single void in an extended solid.] Figure 2 shows the void 
shape changes for the indicated cases. 
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4. EFFECT OF GRAIN BOUNDARY SLIDING ON NONLINEAR STEADY CREEP 


At elevated temperatures, creep of polycrystals involves nonlinear flow within 
grains accompanied by grain boundary sliding which can be modelled by a linearly 
viscous relation; see Ke (1947) , Zener (1948) , and McLean (1957) . The problem of 
estimating the overall creep properties of a polycrystal on the basis of different 
constitutive relations for the grain and the grain boundary has been examined by a 
number of researchers using various models; see, e.g., Zener (1948), Budiansky and 
O'Connell (1976), and Chen and Argon (1979). Recently, Ghahremani (1980a, b) has 
studied a two-dimensional model of creep using a numerical approach. Except for his 
work, other studies do not include the full effect of the essentially periodic 
structure of the grain boundary geometry, and hence the corresponding interaction 
effects. 

In this section we shall examine the creep of polycrystals on the basis of non- 
linear transgranular and linear intergranular creep laws, using a two-dimensional 
(plane) model. 

Figure 3 shows a typical unit cell of dimensions and A^. Within the matrix, 

the flow is governed by constitutive relations (2.3) which, in conjunction with a 

linear creep in bulk, de, , = k do , k = constant, yield 

Kk mm 


^ijki 


—[I 


IjkjL 


n '^ii‘^k£ 1 


^ij'^kj"^ 2 k ^ij^kL’ 


(4.1) 


so that da. . ” C . , 

ij Ijki 

within the grains. 


dcj^j^ holds for the incremental stress, strain-rate relation 
In view of (4.1), Eq. (2,14) must be replaced by 


®ljmn “ ^jk^l^^^Alonn^Z’ 


where, now, N., is 
jk 


(4.2) 
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N ■ — f— ^ -f — ^^6 + — S i. 

jk D L I 2 ® k ^ n+1 

nJ 




2J^ 


(I,. - 6,, 5 

jkmn jk mn 


)} 


•" fc 5j5k>] • 


(4.3) 


where 


1 


r4 _ 2 

‘^k'ki" 

2 

2n J 

- 

2n L2 

^ n+1 ^ 

2J^ 

4. 

1 

ri ?*+ _ _2_ 

/p2r 


2xrnj" 

L2 ^ n+1 

a 




5, a’ a' C, 


To apply Eq. (2.23), we must calculate the quantity Q^(5) for the typical r^^ 
grain boundary segment. For a two-dimensional model, this is easily done and, if 
denotes the center of the segment, and 0 its orientation relative to the x^-axis, 
see Fig. 3, then we obtain 

q (O - e - h( 5 ,e’'), 

. / _ .r, sin sin 

h(±5,6 ) - — ^ , 

7^ (4.5) 

r 

z"" » -—[Si cos 9"^ + ^2 sin s’"] , 
r 

y = “["^1 ^2 

where st^ is the length and t^ the thickness of the r^^ grain boundary segment. Note 
that Eq. (2.21) now becomes 

M 


±CO 


J 's 'r 8ij^ia)l>(5.«'^)h(5.6=)cos{;-(xJ-x=)}dJ*=. (4.6) 

J J s*=l n “0 

P 

T ^ 

Note also that when the thickness t is small relative to the length 1 of a segment, 
then sin y*‘/y^ = 1 in (4.5) . 

Z J H 
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The unit cell shown in Fig. 3 includes a total of 9 grain boundary segments, 
so that M = 9 in Eq. (4.6) and in Eqs. (2.21). For each stress increment (or the 
strain rate increment), we first solve (2.21) to obtain d^ r = 1,2,..., 9. Then 
we calculate the stress increment and update the overall total stress. With this 
stress, we calculate the instantaneous moduli of the matrix from the nonlinear creep 
law (2.3). Equation (2.23) finally yields the overall instantaneous moduli. 

Table 2 shows the geometrical data for the considered unit cell. It is easily 
seen that, in this case, 


= TT [n^^ cos 0 ^ +/3 n^ sin , 

r , to ^ 

y = TT ^ [- sin e + /3 n 2 cos 0 ] , 
where tQ = 3 t/Aj^; note that 


f = 





(4.7) 


(4.8) 


For the numerical calculations, we have assumed 

^ = 1.001 so that K =“ K , 

H 

H = 0 so that T] » 1 , (4.9) 

"22 

n = 3 , and p = — — = 0 (uniaxial tension) . 

‘^ll 

Detailed results are obtained for two cases: (1) tQ = 0.1 which implies that 

f =“ 5.8%. We note that the model considers the linear viscous flow in a rather 
thick band about the grain boundary, and a nonlinear power law with n = 3 (in Eq. 
(2.3)) outside of this band. This model appears reasonable when we observe that 
Instead of the local stress we have used the overall average stress in calculating 
the instantaneous moduli for the grains. 

The results are presented in terms of the following nondimensional quantities: 
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(4.10) 



leads to 

"" ^2222 " "" ^1111 " 2 2 ^ ’ 

(4.12) 

^ ^1122 " ^ ^2211 " 2 ~ 2v ’ 

^ ^1212 ^ ^2121 2v 

In Fig. 4, results are plotted in terms of non-dimensional axial stress and 
strain measures, instead of the effective stress and strain. At stress levels near 
Sq = 1, the lateral strain, C 22 > is positive (extension) and larger than and 

k,Cii 22 is negative for smaller S^. This anomalous result stems from the assumed 
power law creep for the matrix. Another peculiar phenomenon at this stress level 
is that some of the overall moduli are negative ; the shear modulus remains positive. 
Another anomalous behavior for power law constitutive relations has been observed 
by Budiansky £t a^. (1982), in connection with void growth. These authors report 
examples in which, under axial tension larger than the lateral ones, a void in a 
power law matrix is predicted to extend more rapidly laterally than in the axial 
direction. 

The results in Fig. 4 are tentative, as we are now examining this problem in 
more detail. 
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Table 1: Initial and loading conditions 

for considered cases of void growth problems; 
Oq - Bq - 1.0, and fg - .005. 



Table 2i Geometrical data for grain boundary 
configuration in a unit cell. 



133 





































































































134 



135 





A THEORY FOR THERMOVISCOPLASTICITY 
FOR MECHANICAL AND THERMAL LOADING 

E . Krempl 

Rensselaer Polytechnic Institute 
Troy, New York 12181 

and 

E.P. Cernocky 
University of Colorado 
Boulder, Colorado 80302 


Abstract 


A coupled isotropic thermoviscoplasticity theory for small strain 
is proposed. The theory consists of a mechanical constitutive equation and 
a constitutive assumption for the heat equation. These equations are sepa- 
rately postulated but are coupled through their common linear dependence upon 
stress rate and the mechanical strain rate tensors and the time rate of tem- 
perature. The equations depend nonlinearity on the stress and strain tensors 
through the overstress tensor which is the difference between the stress 
tensor and the equilibrium stress tensor (obtained as the loading rate approaches 
zero) and on the absolute temperature. The concept of a yield surface is not 
used and the transition from linear thermoelastic behavior to nonlinear inelas- 
tic behavior is smooth. Extensions of the theory to cyclic loading are under 
development. 

The theory is first applied to conditions of homogeneous deformation 
where the temperature changes in the material are induced by deformation alone. 
For adiabatic conditions numerical experiments (the integration of the coupled 
nonlinear differential equations for the conditions employed in materials test- 
ing using postulated material functions) show that the theory reproduces initial 
thermoelastic behavior (cooling (heating) in uniaxial tension (compression), 
isothermal behavior in torsion) followed by inelastic heating in any state of 
stress during monotonic loading. The amount of deformation induced temperature 
change is negligible unless the loading is very fast. During cyclic plastic 
loading the temperature increase can be considerable and it is shown that the 
predictions of the theory compare very well with experiments performed at room 
temperature on Type 304 Stainless Steel and on a 3.5 Ni-Mo-V steel. 

When large temperature changes are imposed the deformation induced 
temperature changes can be neglected. The numerical experiments involve in this 
case the uniformly changing temperature and the mechanical loading as inputs 
(no heat conduction is allowed). Although other possibilities exist only the 
elastic modulus is assumed to be a function of temperature. The response of 
the model is shown for heating and thermal cycling under mechanical constraint 
(thermal fatigue) and for combined thermal and mechanical cycling of a uniaxial 
bar. It is shown that the response depends on the rate of temperature applica- 
tion and on the temperature at which clamping occurs. 


137 



References 


1) E. P. Cemocky and E. Krempl, Int. J. Solids and Stvuatio’eSj 16, 723-741 
(1980) . 

2) E. P. Cernocky and E. Krempl, J. of Thermal Stresses^ 4^ 69-82 (1981). 

3) E. P. Cernocky and E. Krempl, J. de Meaanique Appliquee^ 5^ 293-321 (1981). 

4) S. L. Adams and E. Krempl, "Thermomechanical Response of 3.5 Nl-Mo-V Steel 

and Type 304 Stainless Steel under Cyclic Uniaxial Inelastic Deformation," 
RPI Report MML-82-5, April 1982. 


138 



